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ABSTRACT 

We report the first measurements of anisotropy in the cosmic microwave background (CMB) radiation 
with the Arcminute Cosmology Bolometer Array Receiver (Acbar). The instrument was installed on 
the 2.1 m Viper telescope at the South Pole in January 2001; the data presented here are the product of 
observations up to and including July 2002. The two deep fields presented here, have had offsets removed 
by subtracting lead and trail observations and cover approximately 24 deg^ of sky selected for low dust 
contrast. These results represent the highest signal to noise observations of CMB anisotropy to date; in 
the deepest 150 GHz band map, we reached an RMS of ^ 8.0 /xK per 5' beam. The 3 degree extent of 
the maps, and small beamsize of the experiment allow the measurement of the CMB anisotropy power 
spectrum over the range £ = 150 — 3000 with resolution of A£ — 150. The contributions of galactic dust 
and radio sources to the observed anisotropy are negligible and are removed in the analysis. The resulting 
power spectrum is found to be consistent with the primary anisotropy expected in a concordance ACDM 
Universe. 

Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

Observations of the Cosmic Microwave Background 
(CMB) provide an unique probe of the Universe at the 
epoch of matter and radiation decoupling. The physics 
of the early Universe can be described in terms of models 
that yield precise predictions for cosmological observables. 
Within the context of these models, observations the CMB 
anisotropy power can be used to place constraints on the 
values of key cosmological parameters (White et al. 1994; 
Hu & White 1996) 

On sub-horizon size scales, gravity driven acoustic os- 
cillations in the primordial plasma give rise to a series 
of harmonic peaks in the CMB angular power spectrum 
(Sunyaev & Zel'dovich 1970; Bond & Efstathiou 1987; Hu 
et al. 1997). A number of experiments have characterized 
the power spectrum up to and including the third acous- 
tic peak (Lee et al. 2001; Netterfield et al. 2002; Halverson 
et al. 2002). These observations have been used to pro- 
duce constraints on cosmological parameters such as the 
total energy density and baryon density with a precision 
of ~ 5 - 10% (Lange et al. 2001; Jaffe et al. 2001; Pryke 
et al. 2002; Abroe et al. 2002). 

On fine angular scales, the CMB power spectrum is ex- 
ponentially damped due to photon diffusion and the finite 
thickness of the surface of last scattering (Silk 1968; Hu 
& White 1997). Observations of the CMB power spec- 
trum in this region can be used to produce independent 



constraints on cosmological parameters. For example, the 
scale of the damping can be used to constrain ^Im and 
flB (White 2001). Recent observations with the Cosmic 
Background Imager (CBI) have provided a first look at the 
damping tail of the CMB and found it to be consistent with 
models motivated by observations on larger angular scales 
(Sievers et al. 2002). On angular scales of a few arcmin- 
utes, deep pointings with the CBI detect power in excess of 
that expected from primary CMB anisotropy (Mason et al. 
2002). In a companion paper to that work, this signal is 
interpreted as being due to the Sunyaev-Zel'dovich Effect 
(SZE) in distant clusters of galaxies (Bond et al. 2002). Al- 
ternative interpretations have been proposed that explain 
the excess power as being due to local structure and non- 
standard inflationary models (Cooray & Melchiorri 2002; 
Grifliths et al. 2002). If the signal is due to the SZE, it's 
precise characterization would provide information about 
the growth of cluster scale structures (Holder et al. 2001; 
Zhang et al. 2002; Komatsu & Seljak 2002). 

The Arcminute Cosmology Bolometer Array Receiver 
(Acbar) is an instrument designed to produce detailed 
images of the CMB in three millimeter-wavelength bands. 
This paper is the first in a series reporting results from the 
Acbar experiment, and describes the CMB angular power 
spectrum determined from the 150 GHz data. In a com- 
panion paper, the Acbar CMB power spectrum and re- 
sults from other experiments are used to place constraints 
on cosmological parameters (Goldstein et al. 2002). Fur- 
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ther papers describing the AcBAR instrument (Runyan 
et al. 2002b), SZE cluster searching (Runyan et al. 2002a), 
and pointed SZE observations (Gomez et al. 2002) are in 
preparation. 

We present a brief overview of the telescope, receiver, 
and site in § 2. Our observing technique, including data 
editing and calibration are described in § 3. The analysis of 
the data is presented in § 4, including several developments 
in the treatment of high sensitivity ground based CMB ob- 
servations. The main results of the paper are presented in 
§5. In § 6, potential sources of foreground emission and 
their treatment in the analysis are discussed. Tests for 
systematic error in the power spectrum are discussed in 
§ 7. In section § 8, we present our conclusions. 

2. ACBAR INSTRUMENT 

The Arcminute Cosmology Bolometer Array Receiver 
(Acbar) is a 16 element 235 mK bolometer array that has 
been used to image the sky in three millimeter-wavelength 
bands. The instrument was designed to make use of 
the Viper telescope at the South Pole to produce multi- 
frequency maps of the CMB with high sensitivity and high 
angular resolution (~ 4 — 5'). 

The Viper telescope is a 2.1m off-axis aplanatic Gre- 
gorian telescope designed specifically for observations of 
CMB anisotropy. A servo controlled chopping tertiary 
mirror is used to modulate the optical signal reaching the 
Acbar receiver. The secondary mirror produces an im- 
age of the primary mirror at the tertiary and therefore the 
motion of the chopping mirror produces minimal changes 
in the illumination pattern on the primary. With this sys- 
tem, it is possible to modulate the 16 Acbar beams 3° 
in azimuth in a fraction of a second without introducing 
excessive modulated telescope emission or vibration. The 
primary is surrounded by an additional 0.5 m skirt that 
reflects primary spillover to the sky. To minimize possi- 
ble ground pickup, a reflective conical ground shield com- 
pletely surrounds the telescope, blocking emission from 
elevations below ~ 20°. 

Between 1998-1999, Viper was equipped with the single- 
element HEMT-based CORONA receiver operating at 
40 GHz. These observations produced a detection of the 
first acoustic peak in the CMB power spectrum (Peterson 
et al. 2000). The Acbar instrument was designed specif- 
ically to take full advantage of the unique capabilities of 
the Viper telescope and was deployed to the South Pole in 
December 2000. 

Acbar makes use of microlithographed "spider-web" 
bolometers developed at JPL as prototypes for the Planck 
satellite mission. The detectors are cooled by a three stage 
closed-cycle ^He-^He sorption refrigerator to 235 mK, at 
which point they are background limited. The fridge and 
focal plane are moimtcd on the 4K cold plate of a liq- 
uid helium and nitrogen cryostat. The signals from the 
bolometers are amplified and sampled with a 16 bit A/D 
at a frequency of 2.4 kHz 

A set of sixteen corrugated feed horns couples the radi- 
ation from the telescope to the detectors. The horns are 
designed to produce nearly Gaussian beams on the sky 
that are all of approximately equal size at all observing 
frequencies. The focal plane array projects onto a 4 x 4 
grid on the sky with a spacing of ~ 15' between array el- 



ements. Behind each beam defining horn there is a filter 
stack that defines the passbands and blocks high frequency 
leaks. The band centers are 150, 220, and 280 GHz with 
associated bandwidths of 30, 30, and 50 GHz, respectively. 
In the 2002 observations, each of the 8 150 GHz channels 
typically achieved a noise equivalent CMB temperature of 
~ 340 /xKa/s. Details of the instrument construction and 
performance arc presented by Rimyan et al. (2002b). 

The South Pole station is located at an altitude of 
~ 2900 m and experiences ambient temperatures ranging 
from — 30°C to — 80°C. In the best winter weather, the per- 
ceptible water vapor has been measured to be ~ 0.2 mm 
(Chamberlin 2001). The high altitude, dry air, and lack 
of diurnal variations result in a transparent and extremely 
stable atmosphere (Lay & Halverson 2000; Peterson et al. 
2002). During the winters of 2001 and 2002, the typical 
atmospheric optical depth for the AcBAR 150 GHz chan- 
nels was measured to be approximately 3%. The entire 
southern celestial hemisphere is available year round al- 
lowing very deep integrations. The combination of these 
unique features and the established infrastructure for re- 
search make the South Pole a nearly ideal site for ground 
based observations of the CMB. 

During an observation, the telescope tracks the position 
of the observed field. Due to the proximity of the telescope 
to the geographic South Pole, the telescope needs to ro- 
tate in azimuth with only small changes in elevation. The 
beams of the array follow a constant velocity 3° triangle 
wave in azimuth with a speed fast enough that the atmo- 
sphere is essentially stationary, and slow enough that the 
beams take several time constants to move across a point 
source on the sky. Subject to these constraints, the chop 
frequency was chosen to be 0.7 Hz in 2001 and 0.3 Hz in 
2002. 

The combination of large chop 3°) and small beam 
sizes (^ 4 — 5') make AcBAR sensitive to a wide range of 
angular scales (150 < £ < 3000), with high £-space res- 
olution (A^ 150). Another unique feature of Acbar 
is its multi-frequency coverage, which has the potential 
to discriminate between sources of signal and foreground 
confusion. The CMB power spectrum we present in this 
work is derived from the 150 GHz channel data. In 2001, 
Acbar was configured with four 150 GHz detectors; this 
number was increased to eight for the 2002 observations. 
An analysis of the 220 GHz and 280 GHz data is underway. 

3. OBSERVATIONS 

To minimize possible pickiip from the modulation of 
telescope sidelobes on the ground shield, we restrict the 
CMB observations to fields with EL > 45°. Fortunately, 
the lowest dust contrast region of the southern sky is cen- 
tered at an elevation of EL~ 55° when observed from the 
South Pole. Several low dust, high declination fields were 
selected for CMB observations. In general, the Acbar ob- 
servations have focused on producing high signal to noise 
maps rather than covering more sky in order to mini- 
mize the sensitivity of the power spectrum to the details 
of the noise estimate. The power spectrum reported in 
this paper is derived from observations of two separate 
fields, which we call CMB2 and CMB5. Each field 
was chosen to include a bright quasar, PMN J0455-4616 
(aj2ooo = 4''55™50.8^ 6j2ooo = -46°15'59") and PMN 
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J0253-5441 {aj2ooo = 2'^53'"29.2^ Sj2ooo = -54°41'51"), 
respectively. As described in §3.3, the pointing model is 
derived from frequent observations of quasars, Galactic 
HII regions and planets. The images of the guiding quasars 
produced during the CMB observations provide stringent, 
independent constraints on beam sizes and pointing accu- 
racy. The two fields are sufficiently separated to be con- 
sidered independent, and the signal correlation between 
them can be safely ignored. 

3.1. The Lead- Main- Trail Scan Strategy 

During the CMB observations, the telescope tracks a 
position on the sky while the chopper sweeps the beams 
back and forth three degrees in azimuth. The motion of 
the chopper introduces systematic offsets into the data 
that must be treated in the analysis. As will be described 
in §3.3, the motion of the chopper causes the beams to 
move slightly 2') in elevation as they sweep in az- 
imuth. This motion modulates the atmospheric emission 
and introduces an ~ 5 mK(T(jir/200 K)(t/0.03) signal at 
an elevation of 45°, where r is the zenith optical depth at 
150 GHz. Although the telescope optics were designed to 
produce a stationary illumination pattern on the primary 
mirror, the illumination on the secondary mirror changes 
substantially with chopper angle. Snow accumulation or 
temperature differences on the secondary mirror will be 
modulated and contribute a second source systematic sig- 
nal that varies with the chopper position. 

The efifects described above are collectively referred to 
as the chopper synchronous offset. Changes in the efi'ec- 
tive temperature of the atmosphere and optics through- 
out an observation can result in a time variation of this 
signal. In order to minimize the chopper synchronous 
offset, each field is observed in a rapid Lead-Main- Trail 
[LMT) sequence where the Main field is led and trailed 
by two identical Lead and Trail observations each with 
half the integration time of the main observation. The 
CMB power spectrum is derived from the differenced map. 
Main- (Lead-|- Trail) /2. If the LMT switching is performed 
faster than the time scale of any drift in the chopper syn- 
chronous offset, then the offset can be eliminated. Com- 
pared to the differencing of two equally weighted fields, 
the LMT observation further removes a linear drift in any 
chopper synchronous offset. In practice, the Lead field is 
tracked for 30 seconds; the telescope is moved -|-3° sec 5 
(in increasing RA) to track the Main field for 60 seconds; 
the telescope is moved another +3° sec 5 to the Trail field 
for 30 seconds of integration. After completing the LMT 
cycle, the elevation of the telescope is decreased by 1' and 
the process is repeated ^ 100 times to build up a con- 
tinuous 2-dimensional map. The CMB2 and CMB5 fields 
actually consist of two largely overlapping sub-fields offset 
by 0.5°sec((5) in RA. The total observation time for each 
field is split roughly equally between the subfields which 
are referred to as CMB2a,b and CMB5a,b. 

3.2. Data Cuts 

In order to minimize the potential contamination of the 
data by systematic errors, we introduced a number of con- 
servative data cuts. All 4 hour observations containing re- 
frigerator temperatures higher than 250 mK cut from the 
data set due to their uncertain calibration. Cryogenics 



fills, refrigerator cycles, computer and telescope mainte- 
nance, and bad weather all limit observation time. In- 
cluding galactic observations and calibrations, the effec- 
tive observing efficiency was ~ 60% for the winter of 2002. 
The effects of cosmic ray interactions with the detectors 
were removed by discarding any chopper sweep contain- 
ing a point varying from the mean of the raw 2.4 kHz data 
samples by > 6(t. Due to the low filling factor of the micro- 
mesh detectors, this results in an insignificant loss of data. 
Before performing the LMT differencing, the data are ex- 
amined for the presence of excessive chopper synchronous 
offsets that indicate the accumulation of snow on the mir- 
rors. The offset level is bimodal between the two extremes 
and the amount of data cut depends only weakly on the 
cut level. Files with snow accumulation are discarded, re- 
sulting in a loss of approximately 40% of the remaining 
data. The details of the cut level and amount of discarded 
data have no significant effect on the resulting maps or 
power spectrum. After implementing the data cuts de- 
scribed above, we retained ^ 400 hours of observation on 
CMB2 during July-August of 2001 and April of 2002. For 
CMB5 we retained ~ 800 hours of data gathered during 
April-July of 2002. As described in §4.1.3, we cut an addi- 
tional ~ 25% of the remaining data that we determined to 
be corrupted by atmosphere induced correlations between 
array elements and adjacent scans. 

3.3. Pointing 

The Viper telescope pointing model is developed from 
observations of bright galactic sources, planets, and 
quasars. Daily observations of the compact galactic 
sources RCW38, RCW57 and MAT6a are used to mon- 
itor pointing offsets and check the pointing model. The 
pointing model contains terms for encoder offsets, colli- 
mation error, telescope flexure, and tilt of the AZ ring. 
When the telescope is pointed at the equator, the detector 
beams sweep out a strip of approximately constant eleva- 
tion under the motion of the chopper. However, at the 
elevation at which the CMB observations ar(^ made, the 
beams sweep out arcs with amplitude of approximately 2 ' 
across the 3° chop. Observations of planets and galactic 
sources at a number of positions in the chopper throw have 
been used to characterize this effect. For each array ele- 
ment, we determine a pointing model that is a function of 
the elevation, azimuth, and chopper encoders. After ap- 
plying the model, we find the RMS pointing error for the 
2001 season to be ~ 1.4' in RA and 30" in declination. 
The majority of the spread in 2001 was due to a malfunc- 
tioning sensor of the chopper angle. After this sensor was 
replaced in January 2002, the pointing RMS was deter- 
mined to be ^ 20" in RA and ^ 20" in declination. Due 
to early concerns about the accuracy of the pointing, we 
chose each CMB field to have a bright ( > 1 Jy at 5 GHz) 
quasar in the middle. In this way, we have an continu- 
ous monitor and independent check of the accuracy of the 
pointing. 

3.4. Beams 

High signal-to-noise maps of planets (Mars and Venus in 
July 2001, and Venus in September 2002) are used to ac- 
curately measure the beam patterns of the array elements. 
The measured maps were deconvolved with a constant 
temperature circular disk corresponding to the planet sizes 
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Fig. 1. — The beams of the 150 GHz array elements determined from observations of the CMB2 guiding quasar in the 2001 season. These 
images represent an average over the entire observation period and include any distortions due to changes in pointing or beamsize. 




Fig. 2. — The beams of the 150 GHz array elements determined from observations of the guiding quasar in the CMB5 field. These images 
represent an average over the entire observation period and include any distortions due to changes in pointing or beamsize. 



(typically small compared to the beams) , to determine the 
instantaneous beam parameters. 

The coadded image of the guide quasar is used to deter- 
mine the effective beamsizes at the center of the map and 
includes any smearing due to drifts in pointing that oc- 
cur over the period in which the data are acquired. Focal 
plane maps for the array in both the 2001 and 2002 con- 
figurations are shown in Figures 1 and 2. In Table 1, we 
list the FWHM beam diameters derived from Gaussian fits 
to the quasar images for both years. These beam param- 
eters were used in the analysis of CMB power spectrum. 
These results are consistent with the beamsizes measured 
on planets and the observed pointing RMS determined 
from repeated observations of galactic sources. Some of 
the beams in CMB5 quasar image appear to be slightly 
elongated in declination; we believe this to be due to the 
accumulation of frost on the dewar window during the 2002 
observations. We estimate the uncertainty in the FWHM 
of the beams to be less than 3%. 

Due to aberrations in the telescope optics, the beam 
sizes increase slightly near the extremes of the chopper 
travel. This effect was characterized by mapping plan- 
ets and galactic sources over a range of pointings offset 
in RA. The telescope beams were simulated with a ray- 
tracing package and the resulting aberrations were found 
to be consistent with the observations. The FWHM of the 
beams typically changes by ~ 5% over the entire chopper 
throw, and the throughputs are found to be conserved. 
The measured beams are found to be very close to sym- 
metric Gaussians and we describe their shapes in terms 
of their small distortions. We parameterize the aberra- 
tions with three parameters: the distortions along RA, 
DEC, and an axis inclined at a 45° position angle. When 
the distortions are small, all three parameters are small 
and can be treated with a perturbation analysis. In sec- 
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tions §4.2 and Appendix C, we describe how changes in the 
beam shape are treated in the computation of the power 
spectrum. 

3.5. Calibration 

The calibration of Acbar is described in detail by Run- 
yan et al. (2002b); here we will present a brief summary. 
The planets Venus and Mars serve as the primary calibra- 
tors for the Acbar CMB observations. Due to the lack 
of a significant atmosphere and dielectric emission by the 
soil, the spectrum of Mars closely approximates that of 
a blackbody. The brightness temperature of Mars varies 
as a function of its distance from the sun with mean He- 
hocentric value of ~ 206 K at 150 GHz. The FLUXES^o 
software package uses the brightness temperature adopted 
by Griffin et al. (1986) and an accurate ephemeris to com- 
pute instantaneous values for the Martian brightness tem- 
perature. Following Griffin & Orton (1993) we assign an 
overall uncertainty to the temperature of Mars of 5%. The 
brightness temperature of Venus is a function of frequency; 
following Weisstein (1996), we fit the brightness tempera- 
ture of Venus with a linear fit to the published data. For 
the Acbar 150 GHz band, we adopt a brightness tempera- 
ture of 297±22K, close to the value of 294 ±22 K reported 
by Ulich (1981) for a band centered about 150 GHz. From 
the South Pole, planets are only intermittently available 
for observation and are always at elevation < 30°. Ob- 
serving them typically requires lowering a section of the 
groundshield and waiting for the sky to clear from the 
zenith to the horizon. For these reasons, it is impractical 
to use planets to regularly monitor the calibration of the 
instrument. 

The HII region RCW38 is bright, compact, and has a 
declination similar to our CMB fields; therefore, it serves 
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Table 1 
150 GHz Beam Sizes 
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channel 


FWHMi(') 


FWHM2(') 


P.A.'*(degree) 


2001 
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2001 
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5.10 


4.81 


176 


2001 


C3 
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4.76 


161 
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2002 


D3 


5.71 


4.90 


68 


2002 


D4 


5.50 


5.45 


121 


2002 


D6 


5.53 


5.25 


7 


2002 


D5 


5.69 


5.15 


87 


2002 


B6 


4.94 


4.40 


111 


2002 


B3 


5.66 


5.20 


69 


2002 


B2 


6.79 


5.34 


105 


2002 


Bl 


5.46 


4.67 


43 



" Position angle from +RA, counterclockwise. 



as an excellent secondary calibrator. The responsivities 
of the AcBAR detectors are a function of the background 
loading and change slightly with elevation. Behind a small 
hole in the tertiary mirror, we have installed a chopped 
thermal load that produces a constant optical signal which 
is used to monitor the detector responsivity as a function 
of elevation. The difference in atmospheric attenuation 
between the planet and RCW38 observations must also 
be taken into account. Sky dips are used to measure the 
atmospheric opacity and correct for the effects of atmo- 
spheric attenuation. Because the emission from RCW38 
is somewhat extended (~ 3' FWHM) and beamsizes of 
AcBAR differ slightly between the 2001 and 2002 obser- 
vations, the integrated flux (to an angular radius of 8') is 
used as our secondary calibration standard. 
The integrated galactic source flux is given by 

c [ T I ^ti I^Gdn /Rp\ /e--c-(«-)\ 

(1) 

where Rp and Rg are the responsivities of of the detec- 
tors when observing the planet and galactic source, Vp 

and Vg are the measured voltages as a function of angular 
position, Tp and tg are the zenith optical depths of the 
atmosphere, 6p and 9g are the elevations of the planet 
and galactic source observations, /p(i^) is intensity of the 
planet, and /{v) is the spectral response of the system as 
measured in the lab with a Fourier transform spectrometer 
as described by Rmiyan et al. (2002b). Both the Mars and 
Venus calibrations produce a consistent integrated flux for 
RCW38 at 150 GHz of Si = 145 ± 7 Jy where the uncer- 
tainty reflects the scatter among the array elements and 
measurements. Including the uncertainty for the flux of 
the planets, the detector responsivities, and integral over 
solid angle, we find Sj = 145 ± 13 Jy which is dominated 
by the uncertainty in the planet brightness temperature. 

We can now use RCW38 to monitor the calibration of 
the instrument and transfer the calibration to our CMB 
observations. Observations of RCW38 are scheduled to 
bracket each ~ 5-hour CMB observation. The zenith op- 
tical depth at 350 /im is measured every 15 minutes with 
a tipping radiometer (Peterson et al. 2002) and is extrap- 



olated to 150 GHz using the observed scaling between the 
measured Acbar 150 GHz opacity and 350 /um tipper op- 
tical depths. After correcting for varying atmospheric at- 
tenuation, the integrated RCW38 flux is found to be con- 
stant over all observations with an RMS scatter of < 4%. 
Because both of the CMB fields are observed at eleva- 
tions similar to RCW38, the responsivity can be consid- 
ered to be constant. Although small, the difference in op- 
tical depth between the RCW38 calibration and the CMB 
field must still be taken into account. The conversion from 
observed bolometer voltage signals to CMB temperature 
differences is given by 

dTcMB ^ Sji / e-TGCscCec) \ 

dVboio /^^£(ilL/(,,)di./T/Gdfi Ve-^^"^''"^^^"^V ' 

(2) 

where B^iy) is the intensity of a 2.73 K CMB blackbody, 
and OcMB is the elevation of the CMB observation. 

Given the stability of the frequent RCW38 observations, 
the statistical accuracy of the calibration transfer from 
RCW38 to the CMB fields is quite good. The uncertainty 
in the calibration is a combination of the uncertainty of the 
planet brightness temperatures, beamsize measurements, 
spectroscopy, and accuracy of the calibration transfer. As- 
suming a common overall uncertainty in the calibration of 
Mars and Venus, we estimate the uncertainty in the tem- 
perature calibration of the observed CMB anisotropy to 
be ±10%. 

4. ANALYSIS 

Algorithms for the analysis of total power CMB data 
are well developed and have been tested on balloon-borne 
experiments (Netterfield et al. 2002; Lee et al. 2001). How- 
ever, the ground-based Acbar experiment is subject to 
constraints which require significant departures from the 
standard analysis algorithms. In this Section, we outline 
the Acbar analysis and highlight its unique features. 

4.1. Maximum Likelihood Map and Noise Weighted 
Coadded Map 
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The first step of the analysis is to produce a map from 
the timestream data. Suppose rf„ {a = l..nt) are rit 
time-ordered measurements of CMB temperature. This 
vector can be separated into the noise component Ua 
and the signal component A^jTj, such that da = 
na + ^j AajTj, or in matrix representation, d = n + AT°. 
Here Tj (j = l..np) are the CMB temperatures on each 
sky pixel j, after being convolved with the "beam" , or the 
experimental response function. For these observations, 
this will be a Lead[-1] Main[-|-2] Trail[-1] three point beam 
function. A is the nt x Hp pointing matrix where A^j = 1 
if a € j (measurement a corresponds to sky pixel j); and 
Aaj = otherwise. If n » AT" (as is the case for the 
AcBAR data), the time stream correlation can be deter- 
mined from the data itself; (ricn/j) = Na/3 ~ {dadp). 

In general, a map T is produced from a linear combina- 
tion of the timestream d; 

T = Ld, 

where L is a rip x coefficient matrix. The map "pixel" 
space correlation matrix is given by 

(T T*) = L(dd*)L* = LNL* + LA(T°T°*)A*L* 

= Cat + Ct. (3) 

4.1.1. Maximum Likelihood Map 

In balloon-borne or satellite experiments, the matrix L 
is usually taken to be 

L = = (A*N-^A)-iA*N-\ (4) 

and the resulting map Tm is the maximum likelihood map. 
Using eq. [3] , the correlation matrix for the maximum like- 
lihood map can be derived: 

(T„T^) = (A*N-iA)-i + (T°T°*). (5) 

It was shown by several authors (Tegmark 1997; Ferreira 

& Jaffe 2000) that a maximum likelihood map is not only 
a good visual representation of the data, but also a lossless 
way to compress CMB information; the power spectra de- 
rived from the map and directly from the timestream will 
have the same uncertainty. 

However, for a ground based experiment like AcBAR, 
the vertical gradient in atmosphere emission restricts the 
observation to constant elevation strips. Furthermore, the 
special location of the Viper telescope (~ 1 km from the 
geographic South Pole) results in maps with little cross- 
linking in declination. In addition, the large angular scale 
structure in the timestream is dominated by atmospheric 
fluctuations. The lack of cross-linking, and the need to 
remove large angular scale (> 1°) sky noise results in nu- 
merical problems in the computation of the maximum like- 
lihood map and an alternative analysis is required. 

4.1.2. Demodulation Analysis 

One can select certain linear combinations of d which 
arc believed to be free from sky noise contamination. In 
this case, the coefficient matrix L is formed by rows of 
sky "modulation" patterns, which are usually taken to be 
edge-tapered sinusoidal functions. This method is known 
as the demodulation analysis, and has been used success- 
fully by several ground based CMB experiments (Netter- 
field et al. 1997; Miller et al. 1999; Peterson et al. 2000). 

Throughout the paper the Greek indices will be used to represent 
"pixels" . 



One of the disadvantages of demodulation analysis is 
that the choice of modulation patterns is arbitrary, and 
could result in the loss of useful information. In other 
words, the choice of L could result in a power spectrum 
with uncertainties considerably larger than the those found 
from the maximum likelihood map (Tegmark 1997). Fur- 
thermore, the modulated data, or the "map" Ld, has lit- 
tle to do with the actual CMB temperature distribution 
on the sky, making it difficult to treat foregrounds and 
test for systematic errors. We performed a demodulation 
analysis of the AcBAR data presented in this paper as an 
early check of our analysis technique. The resulting power 
spectrum was found to be roughly consistent with that 
presented here, although with somewhat larger uncertain- 
ties. 

4.1.3. Noise Weighted Coadded Map 

Due to time varying atmospheric emission, there are 
large signals in the data timestream. These signals are 
easy to remove from the timestream, however, they can- 
not be simply removed from the final map. Therefore, 
it is essential that the corrupted modes corresponding to 
the atmosphere be removed before the data are coadded 
to form a map. We propose to use the cleaned noise- 
weighted coadded map as an intermediate analysis step. 
The operation of the "corrupted mode projection" matrix 
n on the timestream d projects out the undesired modes, 
resulting in the cleaned timestream data d = lid. We 
take n to be a block-diagonal matrix, where the length 
of each block corresponds to a chopper sweep. Each 
block removes certain modes in each data strip. If U 
consists of m columns of linearly independent undesired 
modes, the projection matrix block can be defined as 
Hb = 1 - U(U*U)-iU* (Tegmark 1997). It is easy to 
see IlbU = 0, and lib* = lib- 

With the corrupted modes removed, we can define the 
cleaned noise weighted coadded map as 

Ti = ^ ^ '^adct ■ 

XiUJa is the "weight" associated with measurement a on 
pixel i. The normalization condition requires that 

Aj^Wc = 1. 

The cleaned map can be represented as a linear combina- 
tion of the raw timestream; 

TEEELd, where L = AA*nn. (6) 

Here 17 and A are diagonal matrices. 

In practice, each sample a is weighted on its inverse vari- 
ance after the projection of the corrupted modes, 

timestream sample, and the Latin indices will be used to label spatial 
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To clarify the relation between the maximum likeli- 
hood map and the noise weighted coadded map, we in- 
vestigate the special case where the noise N is diagonal 
= ^ap^ti) ^iid modes arc projected out (11 = I), 

-Li - 2_^J^ia<ia - ^ /I ./n-2 

= [(A*N-^A)-^A*N-M].. 

In this case, T is identical to the lossless maximum likeli- 
hood map. In general, T will not be lossless; however, if 
the off diagonal terms of N are small, T will approach the 
lossless map. This is the case for Acbar. 

From eq. [4], and (d) = (n) -|- AT° = AT°, it is easy 
to show that the maximum likelihood map is an unbiased 
estimate of the CMB sky temperature: (T^) = T°. How- 
ever, this is not the case for T, since 

(t) = KT° ^ T°; where K = LA. (7) 

Therefore the transformation of the theory matrix Ct = 
K(T°T°*)K* in eq. [3] is nontrivial and has to be calcu- 
lated in the coadding process. Explicitly, 

Kij — \i ^ WqIIq,^ (8) 

is the transformation matrix between the sky and our 
"map" space. Hereafter we drop the tilde and simply write 
T as T. The untransformed sky temperature is denoted 
byT". 

Note that the contaminated mode projection matrix II 

is completely arbitrary and can be made to adapt to the 
conditions of a given observation. We take IT to be block- 
diagonal, each block lib removes polynomial functions of 
certain orders in each chopper sweep. The order of the 
polynomial to be removed depends on the atmospheric 
conditions during the observation. In this way, data cor- 
rupted by atmospheric emission on large angular scales 
can still contribute to the determination of the small scale 
power. Quantitatively, the polynomial order is chosen to 
be the lowest such that the residual correlation between 
sweeps in the "stare" (data strip of constant DEC point- 
ing) is less than 5%. If after removing a 10th order poly- 
nomial a residual correlation greater than 5% remains, the 
data are discarded. To avoid bias in noise estimation the 
order of the polynomial to be removed is constant over 20 
minutes duration. The correction matrix K for the map 
is built by coadding the correction matrices for the in- 
dividual scans with the appropriate noise weighting. The 
resulting correction matrix K was used in a series of Monte 
Carlo simulations to show the band-powers themselves are 
unbiased by the adaptive removal of contaminated modes. 

4.2. Noise Matrix 

The precise determination of the noise correlation prop- 
erties of the data is a non-trivial task for high-sensitivity 
ground based experiments. From eqs.[3,6], the calculation 
of noise correlation matrix Cjv for T is straightforward 
once the timestream noise N is known: 

CN{%j} = ^i^3 ^ ^gi^ff y^, A^^^n^^ , (9) 

or in matrix notation. 

Cat = AA*nnNn'nAA. 



This is just a weighted average of IINII*. However, since 
the noise in the timestream is far from stationary, the 
determination of N requires careful treatment (see Ap- 
pendix A). We are confident that the noise correlation for 
the data presented here has been determined to to better 
than 5%. 

4.3. Theory Matrix 

The observed CMB temperature T'{9, 0) is the convolu- 
tion of the real temperature T{9, (j)) with the experimental 
response function B{0, </>). For a LMT subtracted experi- 
ment, B{9, (j)) is the convolution of the beam Bo{0, (j)) with 
the beam-switching pattern 

S(f) = 5{v) - [5{v - fi) + 5{v - r2)]/2, (10) 

where (5's are the Dirac 6- functions. For AcBAR, 
fi = (0,A^) and = (O,-A0) , where = 3° esc 6*. 
The Gaussian temperature anisotropy can be described 
by the angular power spectrum Ci, or the mean square 
of the miiltipole moments: Ci = (a£„ja<?m) • When the 
observations are restricted to small fractions of the celes- 
tial sphere, the spherical harmonic series can be approx- 
imated by two dimensional Fourier series. In this limit, 
the 2-D power spectrum density (PSD) is the product of 
|S(k = £/27r) p and Ce, where B is the Fourier transform of 
B. (Saez et al. 1996). For a given power spectrum and ob- 
servational response function, the theory matrix gives the 
temperature correlation between pixel pairs in the map. 
Since the PSD and the autocorrelation function (ACF) 
form a Fourier transform pair, the theory matrix can be 
calculated efficiently with Fast Fourier transform (FFT) 
technique. 

However, the noise weighted coadded map defined in 
eq.[6] combines sky temperatures measured by channels 
with different beam sizes at different chopper positions. 
Therefore, each pixel has different beam parameters, and 
strictly speaking, the FFT algorithm is not applicable. 
The straightforward c;alculation of the theory matrix and 
its derivatives with respect to band-powers involves a com- 
putationally expensive four dimensional integral for each 
pixel pair. A common way to deal with this problem is 
to average the beam for the entire map, and generate the 
theory matrix from the average beam (Wu et al. 2001). 

We have developed an alternative approach that em- 
ploys a semi-analytic correction formula to treat the vari- 
ation of beam sizes in the map. The beams for each mea- 
surement are coadded accounting for the position of each 
sample and allowing for changes in the beamshape. Dif- 
ferent pixel pairs then receive an individual correction to 
the overall FFT result, depending on the pixel-beam pa- 
rameters. The full derivation of this technique is given 
in Appendix C. Due to the small changes in beamshape 
(~ 5% FWHM), the effect of the beamshape correction 
is small and does not result in significant changes in the 
measured power spectrum. 

4.4. Power Spectrum Estimation 

The first step in the estimation of the power spectrum is 
the determination of a function describing the likelihood 
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that the data are described by a given model. Assuming 
Gaussianity, the likehhood function is 



roc ici 



-1/2 



exp 



(11) 



where C = + C^.. Once C.^^ = LNL*, K = LA, 
and T = Ld arc calculated from the time stream, the 
anisotropy power spectrum can be readily estimated by 
the standard maximum likelihood method. The only dif- 
ference is the replacement of Ct with KCtK*. 

4.4.1. Signal-to- Noise Eigenmodes And Foreground 

Removal 

Following Bimn & White (1997) and Bond et al. (1998), 
the data are further compressed using signal-to-noise 
eigenmode truncation (Karhunen-Loeve transformation). 
With slight modifications, this can be done simultaneously 
with the removal of foreground templates. Our foreground 
removal method is conceptually the same as the constraint 
matrix formalism (Bond et al. 1998; Halverson et al. 2002), 
or the pseudo- inversion of Tegmark (1997). The idea is to 
replace the inverse of Hermitian square root or Cholesky 
decomposition (Cjy^^^) in eq.[A4] of Bond et al. (1998) 
with a non-square "whitening" matrix. This procedure is 
described in detail in Appendix D. The foreground modes 
removed from the CMB2 and CMB5 fields are described 
in § 6 and listed in Tables 3 and 4. 



4.4.3. Decorrelated Band-Powers 

Due to the nature of the experiment (partial sky cov- 
erage, correlated noise in the timestream, etc.), the band- 
power values are not independent. The correlations be- 
tween them are given by the inverse of the Fisher matrix 
(eq. [15]). To get the correct imccrtainty estimate for a 
band-power value we either need to marginalize over all 
the other bands, or perform a linear transformation to 
decorrelate the data points. We chose to decorrelate the 
data points to maintain as much information as possible. 
The decorrelation matrices are not unique: our choice is 
the Hermitian square root of Fisher matrix, because it is 
localized and symmetric (Tegmark 1997). In the new ba- 
sis, the local likelihood function is in canonical form. 

4.4.4. Shape of the Likelihood Function 

It was pointed out by Bond et al. (2000) that an off- 
set log-normal function is a much better approximation to 
the likelihood function than a quadratic one. Failing to 
account for this could lead to serious misinterpretation of 
an experimental result, a most notable example being the 
quadrupole moment measured by COBE. In the decor- 
related basis, the likelihood function near the maximum 
{qb = Qb) can be approximated by an offset log-normal 
function 



lnr(q)=ln£(q)-i^i^^^^e^-^ 



(17) 



4.4.2. Iterative Quadratic Approach to Maximum 
Likelihood Band-Power 

The anisotropy power spectrum is parameterized by the 
band-power q, where 



D, 



27r 



-Ci = ^ qBXBi- 



(12) 



A convenient choice of xbi are "tophat" functions, i.e., 
Xbi = 1 for £ G B: and XBt — tor £ ^ B. The maximum 
likelihood qss arc estimated itcratively using the quadratic 
iteration method (Bond et al. 1998). Near the extrema, 
the log of the likelihood function can be approximated by a 
multidimensional quadratic function. The true maximum 
of the original C can then be found by an iterative process: 



QB 



n+l 



qb' 



p^(f-V)b, 



where 



Vb 



and 



(T*C-1Ct,bC-1T) - Tr(C-^CT,i3), 
Fbb' = 2'^(^ "^CtjbC ^Ct,b')- 



(13) 



(14) 



(15) 



Ct,b is the partial derivative of Ct with respect to qs- 
With eq. [12], it is easy to see 



C' 



t,b 



EdCx 



(16) 



The relaxation factor p is a numerical constant between 
and 1. The converged power spectrum is independent of 
p, however setting p < I prevents the estimated qs from 
"overshooting" and wandering to some local extrema. Em- 
pirically, p = 0.5 worked very well, and the gs's converged 
after a few iterations. 



where the offset log-normal parameters Z are defined as 

Zs =ln(to+a;s). (18) 

Even to the quadratic order, the Fisher matrix is only an 
approximation to the curvature matrix of the likelihood 
function near the extrema. To examine the real shape of 
the likelihood function, we explicitly fit the curvature (Tb 
and log-normal offset xb and compared the results with 
the Fisher matrix (now diagonal) and the offsets suggested 
by Bond et al. (2000) (see their eq. [28]). In general we 
found very good agreement. In a few bands, however, 
the uncertainty indicated by the true likelihood function 
is ~ 15% larger than that derived from the Fisher matrix. 
The uncertainty and log-normal offsets reported in this 
paper are determined from fits to the likelihood function. 

4.4.5. Window Functions 

Our goal is to produce a CMB power spectrum that can 
be used to constrain cosmological models. Window func- 
tions are used to convert a given model power spectrum 
to quantities directly comparable to the band-power mea- 
surements of the experiment. A window function WBe for 
band "B" should have the following property: 



{qb) = Y.^WBtl£)r)i 



where qs is the experimental band-power measurement. 
Knox (1999) derived the appropriate window function for 
the band-power quadratic estimator (Tegmark 1997), 



B' 



This form of window function has been used in several 
other experiments (Halverson 2002; Pryke et al. 2002; 
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Meyers et al. 2002), and should serve as a good approxi- 
mation for maximum likelihood band-power. Making use 

of eqs. [15] [16], it can be shown that the pre-decorrelation 
window functions satisfy the following normalization con- 
dition, 

Y^XBlWBnlli = ^BB' . (19) 

In Section §4.4.3, we described the computation of a set 
of decorrelated band-powers from the Fisher matrix of the 
raw band-powers. This same linear transformation is ap- 
plied to the window functions to produce a set of decorre- 
lated window functions. 

It is not necessary to calculate the window functions for 
all t. However, we would like to sub-divide the qs bands 
into bins with finer I space resolution than the qB them- 
selves. Otherwise in order to satisfy the normalization 
condition eq. [19], the pre-decorrelation window fimction 
reduces to a Kronecker-(5 in band index "i?" , or the tophat 
function given by XBt- This is a reasonable approximation 
for experiments like BOOMERANG where the sky cover- 
age, and therefore intrinsic I resolution, is large enough 
to resolve all the structures in the power spectrum. For 
AcBAR, the I resolution of the experiment is compara- 
ble to the expected structure in and it is essential to 
precisely characterize the dependence of each band power 
on the details of the power spectrum. The decorrelated 
window functions shown in Figure 3 were calculated with 
a resolution in I space of = 30. The first window 
function has significant oscillations with A£ ~ 110. This 
appears to be an result of the close spac;ing of the fields 
in the LMT subtraction. For sufficiently smooth theoreti- 
cal models these oscillations will average to a mean value, 
however, the exact dependence of the first band power 
on arbitrary input models could be complex. Numerical 
tabulations of the window functions are available on the 
AcBAR public web site. 

4.4.6. Monte Carlo Simulations 

In order to test our analysis technique, we performed 
300 Monte Carlo realizations of CMB2 and CMB5 maps 
that are consistent with our observed noise statistics and 
the same underlying ACDM model. These maps were 
processed exactly as the real data, including using the 
same projection of corrupted modes, in order to test the 
ability of our pipeline to accurately determine the input 
CMB power spectrum. A quadratic estimator technique 
(Tegmark 1997) is used to rapidly determine the power 
spectrum for each of 300 realizations. In Figure 4, we show 
the resulting average power spectrum overlaid on the in- 
put model. It is clear that the input model is recovered 
without bias, even in the low £ bins that are sensitive to 
the angular scales on which the mode removal is occurring. 

5. RESULTS 

5.1. Maps 

Using the technique described in § 4.1.3, the LMT sub- 
tracted data are combined to form the noise weighted coad- 
ded maps shown in Figures 5 and 6. Due to computational 
limitations, we choose the map pixelization to be 2.5'. This 

http : //cosmology . berkeley . edu/group/ swlh/ acbar 



is a compromise that allows us to recover the power spec- 
trum with minimal loss of information while keeping the 
time and memory requirements reasonable for a single 16 
processor, 64 GB RAM node of the National Energy Scien- 
tific Computing (NERSC) IBM SP computer. These maps 
have had the two guide quasars and one additional de- 
tected radio source (Pictor A) removed; black pixels mark 
the effected areas. Pictor A is observed to be extended and 
we remove a larger number of pixels around the position of 
this source. These maps have not yet had the undetected 
PMN source catalog and IRAS dust templates projected 
out; however, no other significant sources remain. For the 
purposes of presentation, we have smoothed the pixelized 
map with a 4.5' FWHM Gaussian beam. The lower panels 
in Figures 5 and 6 show the RMS noise in the LMT dif- 
ference maps as a function of position. Due to diflFerences 
in map coverage, the noise varies across the map; in the 
central region, the RMS noise per 5' beam is found to be 
17/xK and for CMB2 and CMB5 respectively. On 
degree angular scales, the S/N in the center of the CMB5 
map approaches 100. 

5.2. Band-Powers 

After projecting out the PMN catalog and dust tem- 
plates, the formalism described in § 4.4 is used to to de- 
termine the angular power spectrum. The resulting band- 
powers are correlated at the ~ 20% level before they are 
decorrelated. The decorrelated band power values are 
plotted in Figure 7 with a fiducial ACDM model; these 
same results are listed in Table 2. The corresponding 
decorrelated band-power window functions are plotted in 
Figure 3. Before decorrelation, the window functions are 
typically anticorrelated to those of adjacent bands at a 
level of ^ 15%. After decorrelation, the anticorrelations 
decrease substantially. The decorrelated band-powers and 
tabulated window functions are available on the Acbar 
experiment public website^^. 

6. FOREGROUNDS 

Sources of foreground emission have the potential to 
contaminate the measured CMB power spectrum. Mul- 
tifrequency observations can be used create templates for 
some foregrounds and remove their effect from the power 
spectrum. At our observation frequencies, thermal emis- 
sion from galactic dust, extragalactic radio sources, dusty 
protogalaxies, and the Sunyaev-Zel'dovich effect are all po- 
tential sources of foreground emission (Tegmark & et al. 
1999). 

The observed CMB fields were chosen to lie in a re- 
gion of minimal dust contrast near the southern galactic 
pole. Templates for the dust emission are determined by 
performing a LMT difference on the IRAS/DIRBE maps 
pubfished by Schlegel et al. (1998). The emission is extrap- 
olated to 150 GHz, using the scaling of Finkbeiner et al. 
(1999). From this analysis, we expect contributions to 
the observed temperature anisotropy power in the differ- 
enced maps of 9/xK2 and 70 /xK^ for CMB2 and CMB5 
respectively. This power is dominated by structure on the 
scale of the map and decreases on smaller angular scales 
(cx The contribution of dust to the observed power is 

therefore expected to be completely negligible across the 
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Table 2 
Joint Likelihood Bandpowers 



t range 


l- 


^+ 








X (/xK^) 


75-300 






187 


6767 


1323 


-45 


300-450 


307 


459 


389 


2874 


605 


-501 


450-600 


462 


602 


536 


2716 


498 


-328 


600-750 


615 


744 


678 


2222 


360 


-391 


750-900 


751 


928 


842 


2300 


355 


-154 


900-750 


921 


1048 


986 


798 


153 


-76 


1050-1200 


1040 


1214 


1128 


1305 


208 


-10 


1200-1350 


1207 


1352 


1279 


583 


130 


36 


1350-1500 


1338 


1513 


1426 


628 


134 


111 


1500-1650 


1510 


1649 


1580 


351 


110 


160 


1650-1800 


1648 


1785 


1716 


248 


99 


149 


1800-1950 


1782 


1953 


1866 


361 


132 


277 


1950-2400 


1946 


2227 


2081 


337 


109 


557 


2400-3000 


2377 


2651 


2507 


307 


275 


1888 



Note. — Decorrelated band-powers g^, uncertainty ob , and log-normal offset xb from the joint likelihood analysis of CMB2 and CMB5. 
The PMN radio point source and IRAS dust foreground templates have been projected out in this analysis. This is the same data shown 
in Figure 7. i± are the multipole £ where the window functions drop to half of their peak values; the effective multipole of the band is 
^eff = X] i{We/i)/ J2(Wi/^)j where the average is done for £- < £ < i^. The window function for the first band is oscillatory, and £eff is taken 
to be the average of the upper and lower £ boundary. 



Table 3 
CMB2 Foreground modes 



RA'' DEC Degrees of Freedom 

PMN J0455-4616 04''55'"51^ -46°15'59" 49 

Pictor A 05''19'"50^ -45°46'44" 81 

SFD dust map*" - - 1 

Other PMN sources - - 102 



""Epoch 2000. 

''Predicted 150GHz emission from DIRBE - calibrated IRAS 100 iJ.m map (Finkbeiner et al. 1999). 

Note. — These modes axe removed from the CMB2 field and do not contribute to the power spectrum. 



Table 4 
CMB5 Foreground modes 



RA DEC Degrees of Freedom 

PMN J0253-5441 02''53'"29'' -54°41'51" 49 
SFD dust map - - 1 

Other PMN sources - - 68 



Note. — 



Same as table 3, for CMB5 field. 
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Fig. 4. — The results of 300 Monte Carlo runs using the measured ACBAR noise correlation. The solid line is the input fiducial ACDM 
model. It is clear that the analysis method accurately reproduces the input power spectrum. 



entire AcBAR angular power spectrum. Projecting out 
dust emission using the IRAS/DIRBE template results in 
the loss of only one degree of freedom and has no signifi- 
cant effect on the resulting power spectrum. 

Extragalactic Radio sources are another potential source 
of foreground confusion. By comparing the positions 
of known radio sources from the PMN 5.0 GHz survey 
(Wright et al. 1994) with the CMB maps, we have de- 
tected 3 radio sources with significance greater than 3a. 
Two of these are the guiding quasars deliberately selected 
to lie in the center of the observed fields. In Tables 3 and 
4, we list the coordinates of these sources and the number 
of pixels removed to clean them from the maps. It is pos- 
sible that many PMN sources lie just below the threshold 
of detection and therefore contribute statistically to the 
observed power spectrum. We created templates which 



include all 170 PMN point sources in the CMB field with 
fluxes greater that 40 mJy at 5.0 GHz and use it to project 
out any possible contribution to the measured power spec- 
trum. In Figure 7, we show the CMB angular power spec- 
trum with and without the projection of the undetected 
point sources and the dust templates. Unless stated oth- 
erwise, all results in this paper are found with an analysis 
that removes these source templates. 

It is also possible that dusty protogalaxies could con- 
tribute to the observed CMB signal. Unlike the radio 
sources, we have no template for the positions of these 
objects. Blain et al. (1998) have estimated the contri- 
bution of dusty protogalaxies to measurements of CMB 
anisotropy and predict a contribution to the 150 GHz 
power spectra of (AT)^ 10 (^/2500)2 /iK^. Therefore, 
dusty protogalaxies are not expected to contribute signif- 
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Fig. 5. — The top panel shows the LMT differenced, atmospheric mode removed, noise weighted, coadded map for the CMB2 field. The 
guide quasar and radio source Pictor-A have been replace with black pixels. The small white circle in the lower left hand corner of the map 
represents the FWHM of the average array element beamsize as determined from the coadded quasar image. The predominance of extended 
structure in the vertical direction is because the extended horizontal structure has been projected out in the atmospheric mode removal. The 
lower panel shows the noise in the LMT difference map as a function of position. 



icant signal to the observed CMB anisotropy. 

The most hkely candidate for contamination of the ob- 
served AcBAR power spectrum is the Sunyaev-Zel'dovich 
(SZE) effect in distant clusters of galaxies. Analytic mod- 
els and numerical simulations generally predict a SZE 
power spectrum with a broad peak at angular scales cor- 
responding to ^ ~ 10^ (White et al. 2002). Deep observa- 
tions with the Cosmic Background Imager (CBI), observ- 
ing in a frequency band of 26 — 36 GHz and at angular 
scales corresponding to ^ = 2000 — 3500, have reported a 
significant detection of excess power (AT)^ = 506li[4g fiK^ 
(Mason et al. 2002). Interpreting this power as the SZE 
leads to a value for the RMS mass fluctuations inside 
spheres of 8 h^^Mpc of as ~ 1.0, close to the upper limits 
allowed by observations of the matter power spectra (Bond 
et al. 2002; Komatsu & Seljak 2002). Similar observations 
with the Berkeley Illinois Maryland Association (BIMA) 



array at a frequency of 28.5 GHz and on angular scales 
closer to the expected peak in the SZE power spectrum 
(4g = 6864) also find excess power, (AT)^ = 202^^^^ /zK^ 
(Dawson et al. 2002) that is marginally consistent with the 
results reported by CBI. If these signals are in fact due to 
the SZE, we can use the spectral dependence of the SZE 
to predict the contribution to the Acbar power spectrum. 
Neglecting relativistic corrections, the observed change in 
CMB temperature in the direction of a galaxy cluster with 
Comptonization y is given by 

Ar^TcMBy(^^^?^-4), (20) 

where x — hv/kTcMB- Therefore, the CMB tempera- 
ture anisotropy power due to the SZE should be a factor 
Ar2(150GHz)/AT2(30GHz) = 0.238 smaller for AcBAR 
than was found by CBI and BIMA. In Figure 7 we show 
the AcBAR and CBI power spectra plotted on top of a 
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Fig. 6. — The same maps as shown in Figure 5, but for the CMB5 field. The signal-to-noise in the center of this map for degree-scale 
fluctuations approaches 100. 
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Fig. 7. — The CMB angular power spectrum with (circle) and without (square) the projection of templates for PMN radio point sources 
and thermal dust. Removing these sources of foreground emission has no significant effect on the measured power spectrum. The bandpower 
values and confidence intervals for the foreground projected power spectrum are listed in Table 2. 



ACDM power spectrum. The marginal detection of excess 
power in the last Acbar band-power is consistent with the 
extrapolation of the CBI excess power to 150 GHz if the 
signal is due to SZE. However, due to the large uncertainty 
on this point, the Acbar results are also consistent with 
the origin of the CBI excess being due to CMB anisotropy 
from local structure and non-standard inflationary models 
Cooray & Melchiorri (2002); Griffiths et al. (2002) or a 
new population of faint radio point sources. 



7. SYSTEMATIC TESTS 

7.1. Gaussianity 

The derivation of the band-power q is based on the as- 
sumption that the signal is Gaussian distributed. This as- 
sumption is tested by examining the distribution of signal- 
to-noise eigenmodes (Appendix D). After the maxi- 
mum likelihood band-powers are found, the total covari- 
ance matrix (theory and noise) in the signal-to-noise basis 
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Fig. 8. — The CMB angular power spectrum derived from CMB2 (square) and CMB5 (circle) field. 



is diagonalized and normalized, such that in the new basis 
the covariance matrix becomes an identity matrix. Then 
the data vector is transformed accordingly and com- 
pared to a Gaussian function with unit variance. When 
the 2133 high signal-to-noise modes are compared to the 
Gaussian distribution, the Kolmogorov-Smirnov statistic 
is determined to he D = 0.0159. The corresponding prob- 
ability for values drawn from a Gaussian distribution to ex- 
ceed this value is PksiD > 0.0159) = 64.7%, and the distri- 
bution of eigenmodes is determined to be consistent with 
Gaussian. This result justifies the use of the band-power 
estimation algorithm, however, we point out that this is 
not a sensitive test for any original non-Gaussianity in the 
map, such as skewness introduced by the SZE or other 
secondary anisotropics. The LMT differencing scheme re- 
duces the skewness, and then the arbitrary eigenvector 
directions in the signal-to-noise eigenmode transformation 
further dilutes any residuals. To test the inflationary pre- 
diction that the primary CMB anisotropy is Gaussian dis- 
tributed and to search for non-Gaussian signatures due to 
secondary anisotropies, a reanalysis of the data is required. 

7.2. Difference Map Power Spectrum 

We performed a number of tests for systematic errors 
in the Acbar CMB power spectrum. These tests include 
"jackknife" tests where the raw data for each field is di- 
vided into two sets that can be subtracted to form a dif- 
ference map. The power spectrum derived from the differ- 
ence map is then examined for significant departures from 
zero that would signal the presence of systematic differ- 
ences between the two halves of the data set. This tech- 
nique provides a powerful method to determine if ground 
pickup, calibration, detector time constants, or other po- 
tential systematic effects contribute significantly to the ob- 
served CMB power spectrum. The difference map band- 
powers at high i are also sensitive to errors in the noise es- 
timate. A positive/negative residual in the highest i band 
would indicate an underestimate/overestimate of experi- 
mental noise. 

Because of the adaptive mode removal scheme described 



in § 4.1.3, the differenced power spectrum is not expected 
to be strictly zero. If Ki and K2 are the correction ma- 
trices defined in eq. [8] for two maps Ti and T2, the dif- 
ferenced map (Ti — T2)/2 will have an expectation value 
of (Ki — K2)T°/2. How completely the maps subtract 
depends on the matrix (Ki — K2). Monte Carlo simula- 
tions are used to produce 300 unique realizations of the 
sky (assuming a ACDM Universe) and noise consistent 
with that measured by the experiment. Difference maps 
are generated for these realizations using the K matrices 
determined from the data. The residual difference power 
spectrum and uncertainty are calculated for the Monte 
Carlo simulations and compared to the observed residu- 
als. Typically, the residual due to mode removal is only 
significant when differencing data taken under very dif- 
ferent atmospheric conditions, and then affects only the 
lowest £ bin. 

Two maps were created by separating the data into two 
halves corresponding to the left-going and right-going por- 
tions of each chopper sweep. The power spectra of the dif- 
ference between these maps should be sensitive to system- 
atic errors due to chopper induced microphonic response or 
detector time constants. Because the left and right going 
sweeps are carried out in exactly the same atmospheric 
conditions and are subject to the same mode removal, 
there is no expected residual in the subtracted maps. The 
power spectrum derived from the L/R differenced data is 
shown in the top panel of Figure 9. The error bars on the 
differenced map power spectra are derived from the Fisher 
matrix. The resulting power spectrum is consistent with 
there being no systematic difference between the maps. 

Next, the data were separated by the azimuth of the ob- 
servations into two halves corresponding to observations 
with the telescope pointing in the 180° azimuth ranges 
toward and away from the nearby Martin A. Pomerantz 
Observatory (MAPO) building. The MAPO building is 
the largest object on the otherwise smooth horizon seen 
by Acbar at the Pole and is expected to be the largest 
source of systematic signal due to the modulation of the 
far sidelobes of the Acbar beams. Therefore, the power 
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Fig. 9. — Power spectra of difference maps produced from three different data jackknifes. The top panel shows the power spectrum produced 
from the difference of maps made with left and right going chopper sweeps. The central panel shows the power spectrum produced from the 
difference of maps made with values of azimuth toward and away from the MAPO building. The bottom panel shows the power spectrum 
produced from the difference of maps made from the first and second halves of the data for each of the observed fields CMB2a, CMB2b, 
CMB5a, and CMB5b. In the central and bottom panels the predicted residuals (cross) and uncertainties are determined by 300 Monte Carlo 
simulations of each difference map and include effects of incomplete subtraction due to atmospheric mode removal. In all cases, the difference 
band-powers (diamond) are consistent with there being no significant signal in the subtracted maps. These tests indicate that the noise 
estimate is accurate and the band-power results are free from significant systematic errors. 



spectrum for the differenced map is a sensitive test for 
ground pickup that is not removed by the LMT differenc- 
ing. The resulting power spectra is plotted in the middle 
panel of Figure 9; the results are consistent with the the 
residual calculated from the Monte Carlo simulation. 

In order to test for systematic changes in pointing, cal- 
ibration or beamsize over the course of the observations, 
we also create a difference map by separating the data in 
time. As described in § 3, each of the CMB fields observed 
actually consists of two subfields offset from one another 
by 0.5° sec (5. Due to changes in beamshape across the 
chopper swing, there are subtle differences between the 
subfield maps that could lead to a non-zero residual sig- 
nal when these maps are subtracted. However, when the 
maps are added, the resulting beams are represented by 
the sum of the beams in the input maps and the power 
spectrum should be unbiased. In order to avoid these com- 
plications, differenced maps are created from the first and 
second halves of the data for each of the four observed sub- 
fields CMB2a, CMB2b, CMB5a, and CMB5b. Because the 
data from which the maps are produced are well separated 



in time and are subject to different atmospheric mode re- 
moval, the predicted residuals for the subtracted maps can 
be significant on large angular scales. The lower panel of 
Figure 9 shows the deference map power spectrum with 
the Monte Carlo determined residuals and uncertainties; 
again, there is no evidence for significant excess power in 
the difference map. 

8. CONCLUSION 

We have used the Acbar receiver to measure the angu- 
lar power spectrum of the CMB at a frequency of 150 GHz 
over the multipole range £ = 100 — 3000. The power spec- 
trum we present is derived from approximately 21 weeks 
of Austral winter observations with AcBAR installed on 
the 2.1m Viper telescope at the South Pole. 

In the course of analyzing this data, we have employed 
new analysis techniques designed specifically for high sen- 
sitivity ground based CMB observations. Monte Carlo 
simulations are used to verify that the analysis method 
accurately recovers the power spectrum without bias. The 
power spectrum we present is robust and has passed strin- 
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Fig. 10. — The CMB angular power spectrum for the CBI-Deep, CBI-Mosaic, and ACBAR experiments. The data are plotted on top of a 
fiducial ACDM model. The shaded green bar shows the expected contribution to the ACBAR power spectrum if the excess power found in 
the CBI-Deep data is due to the SZE. The ACBAR data are consistent with the CBI excess being due to SZE, but provide no signifincant 
new constraints on the source of the signal. 



gent tests for systematic errors. Galactic dust emission 
and radio point sources do not contribute significantly to 
the observed power and are projected out in the final anal- 
ysis. Although dusty protogalaxies cannot be ruled out as 
a source of confusion, the expected contribution to the 
measured power is negligible. Overall, the resulting power 
spectrum appears to be consistent with the damped acous- 
tic oscillations expected in standard cosmological models. 
In a companion paper (Goldstein et al. 2002), the Acbar 
power spectrum is used to place constraints on cosmolog- 
ical models. The power in the highest I Acbar bin is 
consistent with the excess power measured in the Deep 
CBI pointings. At this point, the Acbar data lack the 
sensitivity to place significant constraints on the origin of 
the excess power observed by CBI. 
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APPENDIX 

A. NOISE ESTIMATE 

Each Acbar data file contains observations at pointings in declination. Observations at each declination consist 
of lead, main, and trail pointings in RA which we refer to as stares. Each lead or trail stare has 2ns ("-s left going and 
Us right going) chopper sweeps, while the main stares have Arig sweeps. Each sweep contains samples of the detector 
signals. The goal is to estimate the noise correlation in the M — {L + T)/2 differenced, coadded maps for each file, and 
then use this noise correlation in eq. [9] to derive the final pixel space correlation matrix C^r. In order to get better 
statistics for the noise determination, the noise for each sweep is calculated before the sweeps are coadded so that there 
are as many realizations as possible. 
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A file is divided into x Ug sets of sweeps, each containing (2, 4, 2) sweeps in the (L,M,T) fields on na strips in 
declination. The sweeps in the sets of stares are grouped by their time order so that the first sets contains the first sweeps 

from each stare and the last set contains the last sweeps. The sets of sweeps are combined (using M — (L + T)/2) to 
form a ria,-component data vector for each set of sweeps. The correlations are computed from the LMT subtracted sets 
of sweeps. 

The data for one channel is now reduced to the array dataCna;,/!^,/!^] which we write as Dij,k- Assuming that the data 
are dominated by noise, the following moments can be determined 

m2(i,.i, As) = (A.f.fc X Dj^e^k+As){e,k}, 
m3{i,j,Ad,As) = (A,£,fc x Djj+Ad,k+As) {tM}- 

They are interpreted as correlations within a data strip (mi), among data strips in a constant declination {1712), and among 
different declinations (rns). These are all correlations within the same channel. There are also cross-channel correlations 
for TOi, etc., that arc also largely due to the atmosphere. 

As was pointed out in §4, the pixel space correlation matrix is simply the weighted average of IINn*, or (dd*). One 
could calculate it directly from "cleaned" timestream d, however after the mode removal, the correlation within a sweep 
is not just a function of |i — j | anymore, instead it is a function of both i and j. This prevents the application of Wiener- 
Khinchin theorem (PSD ^ FFT ACF) and would make the algorithm time consuming. Instead, the correlation is 
determined before the projection destroys the "stationary" property of the noise by calculating N = (dd*) first, then 
using (dd*) = n(dd*)n*. In Appendix B, the application of the FFT for piecewise stationary random processes such as 
the AcBAR noise is described. A PSD is by definition positive. Due to limited numbers of realizations the PSD calculated 
from the ACF is not always positive. We therefore parametrize the PSD with the first three terms of the Laurent series, 

P5D = ao + al/-l+a2/-^ 

with the constraint that a, > 0. This expression is found to be an excellent description of the noise PSD. 

In principle, the chopper synchronous offset could introduce correlations between declinations. In practice, these 
correlations are effectively removed by the M — (L + T)/2 differencing. In the vast majority of the data, TO3 is consistent 
with zero. It is also found that when the atmosphere is calm, m2 is also small (< 5%) after LMT subtraction and 
projection of corrupted modes. These cross-sweep correlations are included in the noise correlation matrix, but are not 
found to not have a significant effect on the noise estimate. Nonetheless, approximately ~ 20 — 30% of the data was cut 
because it showed a significant residual correlation (rn2 > 5%) between sweeps after atmospheric removal. 

Although the data used in the analysis are dominated by the instrument noise, atmospheric induced cross-channel 
correlations must be taken into account. Although they are small, the inclusion of the cross channel correlations can lead 
to a noise correlation matrix that is not guaranteed to be positive definite. Including these correlations (and removing 
the negative eigenmodes in the S/N transformation [D]) does not appear to effect the resulting power spectrum, however, 
it does introduce an element of uncertainty to the analysis. To avoid this complication, the band-powers presented in this 
paper arc calculated without including cross-channel correlations. This results in a small uncertainty in the noise estimate 
equal to the largest measured fractional cross-channel correlation. By cutting an additional 10% of the data observed to 
have a correlation exceeding 5%, we can guarantee that the noise estimate is correct to this level. From the distribution of 
channel correlations, we estimate the real error in the noise estimate to be much smaller. The power spectrum presented 
in this paper is insensitive to uncertainty in the noise estimate at the level of a few percent. 

B. CORRELATION AND FFT 

The calculation of the various moments described in Appendix A can be done by Fast Fourier transformation. One 
slight complication is Fourier transformation assumes the input data to be periodic. For a strictly stationary process this 
is not a problem since the correlation usually vanishes after a sufficiently long time lag. For piecewise stationary processes, 
a straightforward application of FFT produces false correlation for lag > n/2 , where n is the number of elements. This 
can be prevented by zero padding. 

Suppose we want to calculate the correlation of two real, n-component stochastic data vectors x and y. When x and 
y are stationary, averaging over indices is ensemble average, 

^ n-A 

dj = {XiVj) = C(A = \i- j\) = _ ^ [XjVj+A + Xj+AVj)- (Bl) 

A faster algorithm involves FFT. To prevent periodicity we pad the two data vectors with zeros up to 2n, and Fourier 
transform; 




18 



Kuo et al. 



The inverse Fourier transform of /*<? is 



FFT-\f*g) - 



1 



(2n)2 
1 



2n 

^ xjykexp 

jku—l 
2n 



- — {k-j) 



■ exp 



n J 



2n 



n 



(2n)2 ^ 

= 2^ ^ Xjyk{Sj+A,k + Sj+A,k±2n + Sj+A,k±4n + ' 



jfe=l 
2n 



2n 



n-A 



The last step uses the fact that yfc = for /c > n. Comparing this with eq. [Bl], wc obtain 



C(A) 



-FFT'\rg + ig*) = 



2n 



■FFT-^[n{rg)]. 



(B2) 



n-A ' ' n-A 

Note that with the correction factor 2n/(n — A), eq. [B2] is exact for < A < n — 1. C(A) were then averaged for 20 
minutes to improve signal-to-noise. 

In the presence of a data cut, for example the guiding quasar, a similar "pad-and-corrcct" procedure can be used, such 
that the data corrupted by the quasar will be ignored in the calculation of correlation function. 

C. THEORY MATRIX FOR NON-UNIFORM BEAMS 

For a given 2-D power spectrum density f*(k), the Fourier transform of the temperature map T[r), or the spectrum is 



5(k) 



P(k) 



B(k)C(k), 



where B is the Fourier transform of the response fimction at f = {x,y), and C a Gaussian random variable: C(k) = 
Cr(k) {(^^) = {C,f) = 1. The realization condition requires that 

P(k)=F(-k); C(k)=C(-k). 

Note that (C*(k)(^(k')) = 2^(k — k'), which is a natural consequence of Gaussian random variables. The realization 
condition does not introduce new correlation because (C*(k)(^(— k)) = (<^*(k)^*(k)) = — (^?(k)) = 0. 

The temperature correlation is 



Ct{12} = (r*(fi)T(f2)) = JJ dkdk'exp{-27rz[-k-fi +k'-r2]}(5i*(k)52(k')) 
= jj dkdk'exp{-27ri[-k-fi +k'-r2]}Bi*(k)B2(k') x ^^^^ISEI(^*(k)C(k')) 



= FFT-^ [^(k)6i*(k)^2(k)](,^,^_,^) . 

What we really want to know for band-power estimation is the derivative of Ct{i2} with respect to the band-power qb- 
This is easily done by replacing P in the equation above with proper shape fimctions determined by XBt defined in eq. [12]. 
For simplicity we made XBe flat in D^. In terms of P the filter functions acquire corrections from the transformation 
between and Ce , hence are not flat in band. We parameterize a Gaussian beam as 



Bo(x,y) cx cxp(--); m; 



a p 



I ( X y 



2\a (3 



In this parametrization, if (a,/?, 7) = (cr, ct, 1), Bq is an axisymmetric Gaussian function with FWHM a-\/81n2. In the 
wavenumber domain, the total experimental response function is 



S(k)=S(k)xexp^-l^ 



2 r^,2 



T{kxOi + kyPY {k^a - kyp) 



+ 



2 2 

Here S(k) is the Fourier transform of LMT switching pattern S(f) defined in eq.[10]. For a pair of pixels (1,2), 
Si*(k)S2(k) = |S(k)|2 exp {-47r2 [{fi, + l)alkl + (m2 + l^ykl + iiza^Gyk^ky] } , 

[(7? + lW + (7l + l)«i] 



where 



Ml 



^J■2 



[(7iVl)/?? + (7l + l)/3i] 
4a2 



1, 
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Here {ax,(7y) are taken to be the average of beam RMS in RA and DEC. To the zeroth order in /x 
Ct{12}{h) ~ Ct{i2}(0) = FFT-^ {p(k)|S(k)|2 exp [-47r\alkl + alkl)]} 



(f=f2-fl) 



Since £ = 4:17^ + ky, if ax = cry, B"^ reduces to the well known beam smearing function = exp(— (7^£^) for a 

symmetric Gaussian beam. 

We can obtain next order approximation of Ct{i2} by taking its partial derivatives in /lx: 



da 



T{12} 



dill 

9Ct{12} 



9Ct{12} 



Ai=0 



jii=0 



(r=r2-f 1) 



ju=0 



= FFT-^ {(-4^V^fc^)F(k)|S(k)p exp [-A^^^alkl + a^^^)] } 
= Fi^T-i|(-47rV^A:2)P(k)|S(k)pexp[-47r2(a2fc2 + ^2^2)]| _ _ 

L ) (r=r2— ri 

= FFT-^ {(-47r2a,a^fc,ft,)P(k)|S(k)|2exp [-4^2(^2^2 ^ ^2^2^] | 



(r=r2-ri) 



etc. 



The first few terms of G 



T{12} 



are 



Ct{12} (m) = Ct{i2} (a* = 0) + ^ 



T{12} 



/it=0/ 



T{12} 




The largest corrections occur at high For I ^ 2000, FWHM ^ 5.5' and a 10% beam distortion (/i = 0.1), the leading 
order correction is on the order of An^a^k^fj, ~ 20% for that specific pixel pair. However, since {ax,<Jy) are the noise 
weighted average beam sizes, the correction goes both ways and the resulting power spectrum does not change significantly. 

D. SIGNAL-TO-NOISE EIGENMODE TRUNCATION AND FOREGROUND REMOVAL 

The constraint matrix formalism is typically used to remove undesired foreground modes (Bond et al. 1998; Halverson 
et al. 2002; Meyers et al. 2002). Combined with signal-to-noise eigenmode analysis (Karhunen-Locvc transformation) 
(Bunn & White 1997; Bond et al. 1998), the foreground removal can be implemented in an alternative algorithm. First we 
construct the foreground mode projection matrices that remove the modes described in Tables 3 and 4. Similar to §4.1.3, 
if U consists of m columns of linearly independent foreground modes, in our pixel space U manifests itself as U = KU 
(eq. [7]). Following Tcgmark (1997), we define the projection matrix P = I - U(U*U)-iU*. 

Next we solve for the eigenvalues uji and eigenvectors for the foreground removed noise correlation C'^r = PCatP. 
The projection of P creates am- dimensional degenerate vector space with eigenvalue zero. Numerically, the eigenvalues 
of the undesired modes are at least 10 orders of magnitude smaller than that of the desired modes. Assume ui = U2 = 
■ ■■ = u)m = 0. The "whitening" matrix is defined as 



W 



i 



-1/2 



-1/2, 



such that W*C'^W = I, where I is the n 



I I 

mi dimensional identity matrix. 



After the transformation of the 



temperature map by the application of W, the foreground modes are "projected out" completely, since eigenvectors 
(ei, 62, ■ ■ ■ e,„) with zero eigenvalues are not included in the construction of W. To proceed with the signal-to-noise 
eigenmode transformation, we construct a fiducial theory matrix that is compared with the noise matrix to decide which 
modes contain useful information (high S/N). The theory matrix Ct = J2b iBdCr/dqE, where the qb's are the band 
powers we wish to estimate. The most reasonable choice of a fiducial power spectrum is a constant flat band-power 
= 10^ /zK^, chosen as a conservative upper bound for D^. The application of W to the fiducial theory matrix and 
temperature map removes the corrupted foreground modes: Ct W*CtW; T ^ T^. = W*T. 
By operating on W*CtW with its orthonormal matrix R, it can be diagonalized; 

W*CtW ^ R*W*CtWR = S = diag(£:fe)- 
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Fig. Dll. — The three panels show the noise weighted coadded map T (top), the foreground removed map Tj, (middle), and the high 
signal-to-noise map T(j (bottom). The top panel is the noise weighted coadded map for the CMB5 field; it shows the bright quasar and the 
artifact generated from mode removal (the blue horizontal stripes adjacent to the quasar). The second panel shows the same map after the 
foreground modes defined in Table 4 are set to zero and ignored in the power spectrum calculation. Other than the bright guiding quasar, 
there are no other clear sources of significant foreground emission. This qualitative statement is supported by the fact that after the removal of 
the guide quasar, the power spectrum derived from this field is unchanged by the removal of the radio source and galactic dust templates. The 
bottom panel shows foreground removed map constructed from only the high S/N eigenmodes. Clearly, the S/N eigenmodc transformation 
preserves the basic features of the map. 



In this new basis, the theory matrix consists of diagonal elements of signal-to-noise values Ek- Modes corresponding to 
small £k do not significantly contribute to the likelihood function and only the Ug high S/N modes are retained. For 
the CMB2 field the number of modes was decreased from 9000 to 2500, and for CMB5 the reduction was from 5000 to 
2000. Eliminating these "noisy" modes greatly reduced the computational resources required by the quadratic iteration 
described in §4.4.2. The first columns of R form a mi x Ug transformation matrix Ri. The foreground subtracted, 
high S/N mode data are = R|W*T; and the noise correlation matrix is the dimensional identity matrix. The 
resulting CMB power spectrum is unchanged for values of the cut-off S/N ranging from = 0.1 to 0.01; this is a direct 
consequence of the conservative choice of E^. The results presented in this paper were calculated using £c = 0.05. 

To demonstrate the foreground removal and Karhunen-Loeve transformation in effect, it is useful to reconstruct the 
pixel-space map from the lower dimensional and data vectors, 

T; = W0T^; 

- W0RiT„ 

where = diag(a;/t), (fc — m+1, , m + 2, . . . , n). T'^ is the foreground-removed map, and T's is the foreground-removed 
map reconstructed only from high signal-to-noise modes. TV and T's are both n dimensional, can be directly compared 
with the raw coadded map T. The foreground removed maps for the CMB2 and CMB5 fields are shown in Figures 5 and 
6. For comparison, the raw coadded map, the foreground-free map, and the high S/N eigenmode maps for the CMB5 
field are shown in Figure Dll. 
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